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INTRODUCTION 


Ionospheric  Time-Delay 
Algorithm  for  Single- 
Frequency  GPS  Users 

JOHN  A.  KI  OHK  HAK,  Senior  Member.  IITU 
Aii  I  oreo  (icopliysics  Lalniralory 


The  goal  in  <lcsi}>ui»K  an  ionospheric  time-delay  correction 
algorithm  for  Ihe  siimle-frtqueney  global  positioning  system  user 
w as  to  include  the  main  realurvs  of  the  complex  lielias  im  ol'  the 
ionosphere,  yet  require  a  niininium  of  coefficients  ami  user 
computational  time,  while  still  yielding  an  rms  correction  of  at  least 
50  percent.  The  algorithm  designed  for  this  purpose,  and 
implemented  in  Ihe  <>I’S  satellites,  requires  only  eight  coefficients 
sent  as  part  of  the  satellite  message,  contains  numerous 
approximations  designed  to  reduce  user  computational 
requirements,  yet  preserves  the  essential  elements  required  to  obtain 
group  delay  values  along  multiple  satellite  viewing  directions. 


AoIIio'n  .kljksv  loil'kplk'lk  INlSsks  I  >t\  IMtol.  Nil  I  tore  (  k-ophv  'lev 
i.iil'iH.itoiv  .  ll.Okeom  Ao  lioco  like.  Hcjlnkl.  M  \  stlllll 
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The  U.S.  Department  of  Defense  has  developed  an 
advanced  satellite  navigation  system,  called  the 
NAVSTAR  global  positioning  system  (GPS),  which  will 
allow  users  to  measure  range  and  range  rate  information 
simultaneously  from  four  satellites  to  determine  user 
position  and  velocity  |l  5|.  The  primary  satellite  system 
operating  trequency  is  1.575  GHz.  liven  at  this  relatively 
high  frequency,  the  Earth’s  ionosphere  can  retard  radio 
waves  from  their  velocity  in  free  space  by  nto'c  than 
300  ns  on  a  worst  ease  basis,  corresponding  to  range 
errors  of  l(K)  m.  Liven  monthly  average  diurnal  maximum 
ionospheric  time-delay  values  can  he  as  high  as  150  ns. 
or  50  m  of  range  error,  when  a  GPS  satellite  is  observed 
at  low  elevation  angles. 

In  order  to  nearly  eliminate  ‘lie  effects  of  the  liartlTs 
ionosphere  on  users  of  the  GPS.  a  secondary  frequency. 
1.227  GHz.  was  specifically  incorporated  into  the  system 
to  allow  users  to  automatically  correct  for  the  effects  of 
both  the  range  and  range  rate  errors  induced  by  the 
ionosphere.  Some  users  of  the  GPS.  however,  may  not 
have  a  requirement  for  nearly  complete,  automatic 
correction  lor  ionospheric  range  and  range  rate  errors  as 
the  two  frequency  GPS  provides,  l  or  these  users  an 
algorithm  was  designed,  using  eight  coefficients, 
transmitted  as  part  of  the  satellite  message,  to  provide  a 
correction  for  approximately  50  percent  rms  of  the 
ionospheric  range  error.  Corrections  for  ionospheric  range 
rate  errors  for  a  single  frequency  user  are  not  practical  by 
modeling  techniques,  due  to  the  impossibility  of 
predicting,  except  in  a  statistical  manner,  the  small 
undulations  in  the  ionosphere  which  produce  range  rate 
errors  on  time  scales  of  a  few  seconds  to  minutes.  The 
goal  of  a  50  percent  nils  collection  for  the  ionospheric 
algorithm  was  arrived  at  somewhat  arbitrarily  as  a 
compromise  between  number  of  coefficients  required  io 
be  sent  as  part  of  the  satellite  message  and  the  realization 
that  even  a  state  of  the  art  ionospheric  model,  requiring 
many  coefficients,  would  provide  only  a  70  to  SO  percent 
nils  correction  to  the  ionospheric  time  delay.  In  some 
comparisons  of  the  algorithm  against  actual  ionospheric 
total  electron  content  (i'liO  data  presented  here,  and 
additional  comparisons  given  by  Leess  and  Stephens  |b|. 
it  is  shown  that  the  goal  of  a  50  percent  rms  correct  ion 
for  ionospheric  time  delay  by  use  of  the  algorithm 
described  here  has  been  met. 

CHARACTERISTICS  Of  IONOSPHERIC  TIME  1)11  AY 

The  important  parameter  responsible  for  ionospheric 
time  delay  is  the  total  number  of  electrons  encountered  In 
the  radio  wave  on  its  path  from  each  satellite  to  the  GPS 
system  user.  This  TliC  is  a  function  of  many  variables, 
including  long  and  short  term  changes  in  solar  ioiii/inc 
ilux.  magnetic  activity,  season,  time  of  day.  user  location 
and  viewing  direction.  An  example  of  the  geographic 
variation  of  monthly  average  vertical  time  delay  at  the 
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primary  Ci  PS  system  operating  frequency  IA  is  shown  in 
l  isi.  I  lot  20(H)  U  I .  At  this  time  of  day  the  maximum 
values  of  ionospheric  lime  delay  occur  in  the  American 
longitude  sector.  As  illustrated  in  l-'ig.  I.  the  world  wide 
average  maximum  values  of  TKC  do  not  occur  at  the 
geographic  equator,  but  at  latitudes  approximately  plus 
and  minus  15°  on  either  side  of  the  geomagnetic  equator. 
Note  that  the  time  delay  contours  plotted  in  Pig.  I  are  for 
a  satellite  at  the  zenith.  At  the  minimum  design  elevation 
angle  for  the  GPS  system  of  5‘  the  ionospheric  time  delay 
values  will  be  approximately  three  times  higher  than 
those  shown  in  Pig.  I . 

Many  studies  have  been  made  of  the  behavior  of 
I IX'.  and  the  characteristics  of  this  ionospheric  parameter 
are  now  fairly  well  known,  especially  in  the  midlatitude 
regions  of  the  world,  where  most  single-frequency  GPS 
users  are  likely  to  be  located.  The  behavior  of  TKC'  in  the 
equatorial  and  high  latitudes  is  less  well  known. 

Basically,  the  TKC  behavior  in  many  parts  of  the  world 
has  a  diurnal  maximum  near  1400  hr  local  time,  with  the 
standard  deviation  of  the  TKC  about  the  monthly  average 
value  for  any  given  daytime  hour  generally  between  20 
and  2.s  percent  of  the  mean  value.  The  midlatitude 
ionosphere  generally  has  smoother  latitudinal  gradients 
than  either  the  low  or  high  latitude  regions;  thus,  a 
single-frequency  GPS  user,  in  the  midlatitudes,  has  the 
additional  advantage  of  both  a  smoother  and  more  well 
studied  ionosphere! 

With  these  general  facts  in  mind,  an  algorithm  was 
designed  to  give  best  fit  to  the  large  daytime  values  of 
monthly  average  1 IX'.  and  to  accept  any  difference  from 
monthly  average  TKC  behav  ior  as  part  of  the  residual 
error  lot  a  Gi’S  sx-tem  single-frequency  user.  Since  the 


correlation  distance  |7|  and  correlation  time  |X|  of  TKC 
deviations  from  monthly  average  conditions  are  small, 
and  each  GPS  satellite  message  is  normally  updated  no 
more  often  than  once  each  day.  no  attempt  was  made  to 
make  corrections  for  short  term,  of  the  order  of  a  few 
hours,  variability  of  the  ionospheric  I  PX'.  Also,  since  the 
correlation  of  TKC  daily  values  with  day-to-day  solar  flux 
values  is  low.  no  attempt  was  made  to  tijnlate  the 
algorithm  on  a  daily  basis. 

FORM  OF  THE  IONOSPHERIC  TIME-DELAY 
ALGORITHM 

The  algorithm  developed  is  a  compromise  of  several 
factors.  These  include  I)  user  computational  complexity, 
2)  present  state  of  knowledge  of  temporal,  diurnal,  and 
geographic  variations  of  TKC.  3)  number  of  coefficients 
available  in  the  satellite  message  for  the  ionospheric 
correction  algorithm,  and  4)  likely  single-frequency 
system  user  geographic  operational  areas. 

The  algorithm,  incorporated  in  the  GPS  system  for 
single-frequency  users,  consists  of  a  cosine  representation 
of  the  diurnal  curve,  allowed  to  vary  in  amplitude  and  in 
period,  with  user  latitude.  Fig.  2  shows  how  a  simple, 
positive  half  cosine -shaped  curve  has  been  made  to  fit  a 
typical  monthly  average  TKC  diurnal  variation  from  the 
station  on  Jamaica. 

Pig.  2  illustrates  the  lour  potential  parameters  to  this 
simple  cosine  form  of  time-delay  algorithm  for  the  single- 
frequency  GPS  user.  They  are  I)  the  night-time  constant, 
or  de,  term,  2)  the  amplitude  of  the  cosine  term,  3)  the 
phase  of  the  cosine  term,  and  4)  the  |wriod  of  the  cosine 
term.  Note  that  the  period  of  the  rms  best  fit  is  not  24  hr, 
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LOCAL  TIME 

I'ii;.  2  lA.impIo  of  actual  monthly  average  lime-ilclay  data,  along 
with  cosine  model  fit  to  data. 


but  is  whatever  period  is  required  to  minimize  the  rms 
difference  between  the  actual  monthly  mean  diurnal  curve 
and  the  cosine  tit.  This  equivalent  period  is  generally 
significantly  longer  than  24  hr. 

A  study  of  the  Iwhavior  of  the  do  and  phase  over  a 
wide  range  of  latitudes  showed  that  both  these  terms 
could  he  made  constant  with  little  overall  algorithm  rms 
error  increase.  Thus,  the  de  term  was  made  a  constant 
5  ns  and  the  phase  was  set  at  a  constant  1400  hr  local 
time. 

The  remaining  two  parameters,  the  amplitude  of  the 
cosine  term  and  its  period,  are  functions  of  latitude, 
ordered  by  geomagnetic,  rather  than  geographic,  latitude, 
and  are  represented  in  the  algorithm  by  third  degree 
polynomials.  The  coefficients  of  these  polynomials  are 
transmitted  as  part  of  the  satellite  message.  The  algorithm 
coefficients  were  computed  from  an  empirical  model  of 
world-wide  ionospheric  behavior  derived  by  Bent  |9|,  for 
each  10  day  period  of  the  year  and  for  several  values  of 
average  solar  flux  conditions.  Coefficients  transmitted  by 
the  GPS  satellites  are  updated  once  each  ten  days,  or 
sometimes  more  frequently,  if  the  live  day  running  mean 
solar  flux  changes  by  a  large  amount  during  that  period 
jpf  time. 

•  Simple  cubic  fits  to  both  the  amplitude  and  the  period 
terms  are  not  sufficient  to  adequately  represent  the  actual 
gradients  in  the  equatorial  anomaly  region.  However,  the 
variability  <»i'  equatorial  anomaly  location  and  its 
magnitude  are  large,  not  well  represented  in  the  original 
model  from  which  the  ali  orithm  coefficients  were 
derived,  and.  few  actual  I  PX'  measurements  were 
available  in  this  latitude  region,  when  the  Bent  model 
was  constructed.  Therefore,  no  attempts  were  made  to 
use  an  a  priori  shape  in  the  algorithm  for  the  TKC 
behavior  in  the  low  latitude  ionosphere. 


GEOMETRIC  APPROXIMATIONS -USED  IN  THE 
IONOSPHERIC  TIME-DELAY  ALGORITHM 

The  actual  cosine  representation  of  ionospheric  time 
delay  is  much  simpler  than  the  geometrical  calculations 
required  to  find  the  appropriate  geographic  location,  time, 
and  zenith  angle  for  the  required  time-delay  algorithm 
calculation.  Since  the  GPS  satellites  are  nearly  always 
viewed  at  oblique  elevation  angles,  the  overhead 
ionosphere  is  not  the  appropriate  one  for  a  single- 
frequency  observer  to  use.  The  THC’  must  be  found  at  the 
geographic  point where  the  ray  path  from  each  satellite 
intersects  the  mean  ionospheric  height,  rather  than  at  the 
user  location.  This  point  is  taken  at  a  mean  vertical 
height  of  450  km  along  the  path  to  each  satellite.  T  he 
obliquity,  or  slant  factor,  also  must  be  calculated  for  the 
mean  ionospheric  location,  and  the  vertical  T  EX'  obtained 
from  the  algorithm  at  the  mean  ionospheric  location  must 
be  multiplied  by  ibis  slant  factor.  Finally,  since  the  ITX' 
is  best  ordered  in  geomagnetic,  rather  than  geodetic 
coordinates,  a  conversion  to  geomagnetic  latitude  also  is 
required.  If  the  exact  form  of  the  calculations  for  all  the 
necessary  geometry  were  to  be  used,  the  computational 
assets  would  be  excessive.  Therefore,  in  the  derivation  of 
these  geometrical  functions,  simplifying  assumptions  were 
made  in  many  eases  to  reduce  the  calculation  complexity. 
A  discussion  of  each  of  the  steps  required  in  the 
geometrical  calculations,  and  the  simplifying  assumptions 
made,  follows. 

Earth  Angle 

T  he  TEXT  is  computed  at  the  point  where  the  rav  from 
the  satellite  intersects  the  ionosphere,  rather  than  at  the 
observer.  For  an  observer  located  near  the  earth's  surface, 
looking  at  a  satellite  at  5°  elevation,  the  ionospheric 
intersection  location  is  approximately  14°  of  earth-entered 
angle  from  the  observer.  Since  the  gradient  in  FIX'  can 
be  large  over  Fartb  angles  much  less  than  14"  we  must 
find  the  actual  point  where  the  satellite  intersects  the 
ionosphere,  rather  than  simply  using  the  observer's 
location.  T  he  exact  Faith  angle  representation,  as  well  as 
the  approximate  form,  are  shown  in  Fig.  .1.  versus 
satellite  elevation  angle.  The  approximate  form  is  less 
than  0.2'  in  error  for  all  elevation  angles  greater  than  10  . 
and  is  only  0.4°  and  0.  V1  Emm  the  exact  form  at  5°  and 
0"  elevation,  respectively.  Now  that  we  have  a  useful 
form  of  the  Faith  angle  we  can  compute  the  coordinates 
of  the  ionospheric  intersection. 

Computing  tho  Moan  Ionospheric  Location 

The  actual  form  of  the  equations  to  obtain  the 
ionospheric  intersection  coordinates,  as  well  as  the  Hat  - 
Faith  approximation  used,  are  shown  in  Fig.  4. 
Fortunately,  at  the  low  and  mid-latitudes,  any  errors  in 
using  the  approximate  form  to  calculate  the  ionospheric 
intersection  coordinates  is  small.  To  keep  away  from  the 
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largo  errors  iliis  approximation  would  yield  at  the 
extremely  high  latitudes,  if  the  ionospheric  latitude  is 
calculated  to  he  greater  than  75°.  the  latitude  is  defined  as 
75".  This  approximation  is  reasonable  because  little 
Mngle-lrei|uency  user  activity  is  expected  above  75“ 
latitude,  and  the  TliC  is  generally  low  at  those  very  high 
latitudes  The  few  actual  TIT'  measurements  available  in 
the  polar  cap  region  show  high  variability,  which  is  not 
well  represented  by  currently  available  ionospheric 
models. 


Conversion  from  Geodetic  to  Geomagnetic 
Latitude 

Ionospheric  time  delay  is  best  ordered  by 
geomagnetic,  rather  than  geographic,  latitude  Thus,  it  is 
necessary  to  perform  a  transformation  from  geodetic  to 


geomagnetic  latitude.  This  transformation  from  geodetic 
to  geomagnetic  latitude,  assuming  that  the  earth's 
magnetic  field  can  be  represented  by  an  harth  centered 
dipole,  is  given  by 

sin  <l>  -  sin  (!>  sin  <!>/>  1 II  cos  <|>  cos  «!>/»  costX  X/>) 
where 

«!>/>  -  7X  5  N 
X/>  -  2lM  .(T  I  .. 

The  approximation 

(|>  —  tf)  4  1 1 .6"  cos  (X  —  201) 

represents  the  exact  form  to  within  1“  at  all  geomagnetic 
latitudes  cquatorward  ol  40°.  and  is  within  2“  up  to  h5“ 
either  side  of  the  geomagnetic  equator.  Over  the  entire 
CONUS  region  the  error  in  the  approximate  form  is  less 
than  I". 


finding  Local  Time 

Given  the  approximate  universal,  or  Greenwich,  time 
<UT)  and  the  approximate  longitude  of  the  ionospheric 
point,  the  local  time  at  the  ionospheric  point  is  simply:  / 
=  X,/!5  +  U  l  .  If  r  comes  out  greater  than  24  hr.  simply 
subtract  24  hr  to  keep  the  resultant  time  between  0  and 
24  hr. 


I  he  Obliquity  fac  tor 

The  vertical  time  delay  at  the  suhionospheric  point 
must  he  multiplied  by  an  obliquity,  or  slant,  factor 
defined  as  the  secant  of  the  zenith  angle  at  the  mean 
ionospheric  height.  An  average  ionospheric  height  of 
550  km  is  assumed,  fig.  5  shows  both  the  exact  form  of 
the  slant  factor  calculation  as  well  as  the  approximate 
form  used  in  the  algorithm.  This  approximate  form  is 
within  2  percent  of  the  exact  value  for  all  elevation 
angles  greater  than  5°. 
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7)  Then  compute  the  ionospheric  time  dcluv: 

7,ovo  -  /  x  Is  x  K)  11  !  v  »!»;;, 

|  n-  U 

■('-T-S)' 

where 

2it</  -  SO  UK)) 

A  -  - V -  • 

V  «  «|>» 

Ud  Hm  im 

-  _  II  -  l) —  -  -  -  -  -  ■  ------- 

Note:  7|()N()  is  referred  to  the  L I  frequency.  It  the  user 
requires  the  ionospheric  time-delay  correction  on  the  1.2 
frequency,  the  correction  term  must  be  multiplied  hv  the 
constant  1.65. 


SUMMARY  OF  ALGORITHM  EQUATIONS 

T  he  form  of  the  ionospheric  time-delay  correction 
algorithm  as  given  in  the  GPS  user  "Interface  Control 
Document”  1 10|  is  summarized  in  the  following 
equations.  In  the  actual  implementation  of  the  ionospheric 
correction  algorithm  the  units  used  in  angular  measure  are 
semicircles,  and  time  is  expressed  in  seconds. 

Given.  User  approximate  geodetic  latitude  ‘l>„, 
longitude  X„.  elevation  angle  A\  and  azimuth  A.  to  the 
particular  OPS  satellite  for  which  you  wish  to  calculate 
the  ionospheric  time  delay.  Also  given  are  the 
coefficients.  o„  and  ft,,  which  are  transmitted  as  part  of 
the  satellite  message. 

Note.  All  angles  are  in  units  of  semi-circle;  time  is 
in  seconds. 

1 )  Calculate  the  Liarth-centered  angle 

0.01.17  ... 

i|)  ^ - - -  0.022  (semicircles). 

/•:  ■)•  o.  1 1 

2)  Compute  the  subionospheric  iatitude: 

<l»,  ~  <!>„  )■  i|t  cos  A  . 

If  <l>,  >  t  0.416.  then  <1*,  -  t  0.4 1 0,  If  <l>,  < 

-  0.416,  then  ‘I*,  —  -0.416. 

3)  Compute  the  subionospheric  longitude: 


EXAMPLE  CALCULATION 


As  an  example,  assume  a  station  located  at  40.0“ 
North.  1(K).0°  West,  viewing  a  GPS  satellite  at  20° 
elevation  and  210°  azimuth.  The  ionospheric  algorithm 
coefficients  transmitted  at  this  time  from  the  CiPS 
satellites  are  u„  -  5.82  *  10  s,  1.49  X  10  \  1.79 

x  10  .  Oandp.,  =  1.43  x  10'.  0.  -3.28  x  10'. 
1.13  x  10'.  The  time  is  2045  UT. 

T  he  results  for  this  example  are 

1 )  i|i  —  0,03996  semicircles  (7.2°) 

2)  <l>,  =  0.215  semicircles  (38.7"  N) 

3)  X,  ---  -0.6399  semicircles  <124.795°  \V) 

4)  ll>„,  -  0.2509  semicircles  (45, 16“  N) 


5)  i  =  50.7(H)  s 

6) /  =  2  176 

7)  /'|ON<)  =  77.6  us  (23.3  m). 

'I'his  example  was  chosen  to  represent  an  actual 
coefficient  set  during  a  spring  equinox,  high  solar  flux 
period. 


SUMMARY  Ol  ALGORITHM  CHARACTERISTICS 


4)  Lind  the  geomagnetic  latitude: 

<|i„,  --  <l>,  t  0.064  cos (X/  -  1.617) 

5)  land  the  local  lime: 

t  -  1.32  x  |();  CiPS  time  (see). 

If  I  >  8(>  400.  use  l  -  I  86  l(K).  If  (  -■  O.  add 
86  400. 

6)  Compute  the  slant  factor 

/  -  1.0  +  16.0  x  (0.53  -  /.)'. 


The  ionospheric  time-delay  algorithm  for  the  single- 
frequency  GPS  user  was  designed  to  minimize  user 
computational  complexity,  yet  still  retain  the  goal  of  a  50 
percent  mis  ionospheric  range  error  correction  A  major 
part  of  the  reduction  in  user  computational  complexity 
was  accomplished  by  using  approximations  to  the 
geometrical  calculations  involved  to  reduce  the  number  of 
trigonometric  calculations  which  can  he  computationally 
time  consuming.  In  fact,  the  algorithm  itself,  the  simple 
cosine  representation  of  a  diurnal  curve,  was  truncated  to 
the  first  two  terms  of  the  cosine  expansion  in  order  to 
save  computational  time,  and  incidentally,  slightly 
improve  the  algorithm  fit  to  the  actual  data. 
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INITIAL  TESTS  OF  IONOSPHERIC  ALGORIIHM 
ACCURACY 

Tlu*  ionospheric  time -delay  algorithm  was  tested 
against  actual  THC  from  IS  different  stations  lot  a  total 
of  WU  station-months  of  data  covering  a  representative 
|>ortion  of  the  globe  and  for  both  solar  maximum  and 
minimum  conditions,  big.  6  shows  the  distribution  of 
stations  and  times  from  which  data  were  used  for  model 
tests.  Notice  that  a  significant  portion  of  the  data  is  from 
regions  other  than  the  midlatitudes,  and  that  mans  parts 
of  the  world  are  represented  in  the  available  data.  Also, 
much  of  the  data  Was  from  a  period  of  high  solar  activity 
when  absolute  values  of  THC  are  high. 
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l;ij!  (>  Distribution  ol  limes  liiuI  UkuiImin  when  montlih  axcru^c  TI:C 
data  were  a  vail  aide  tor  algorithm  leMing. 


Summary  statistics  of  ionospheric  time  delay  for  each 
station  for  all  available  hours  and  seasons  were  computed 
in  units  of  mis  error  in  nanoseconds  at  1 .6  Gil/.  Also, 
the  residual  mis  error  for  each  station  was  computed  after 
correction  by  use  of  the  time  delay  algorithm.  These 
summary  results  are  shown  in  big.  7,  with  the  dashed 
line  indicating  the  uncorrected  ionospheric  time  delay, 
and  the  solid  line  indicating  the  residual  time  delay  after 
correction.  T  wo  stations  arc  shown  in  more  than  one 
column  in  big.  7  because  the  actual  ionospheric  TLX'  data 
available  from  those  stations  came  from  more  than  one 
satellite,  and  from  different  time  periods. 

T  he  comparison  of  the  algorithm  against  actual 
monthly  mean  T  HC  data,  for  the  entire  490  station 
months  of  data,  gave  an  overall  rms  residual  error  of  78 
percent.  In  order  to  find  the  rms  error  of  the  model 
against  actual  individual  data,  rather  than  against  monthly 
aserage  ITX'  data,  the  78  percent  algorithm  error  in 
representing  (he  monthly  mean  THC.  was  combined  in  a 
root  of  the  sum  ot  the  square  error  (rss)  sense  with  the  25 
percent  rms  deviation  of  individual  data  from  the  monthly 
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big  7  Kosulis  ol  lomisplierie  lime  delay  algorithm  by  station. 


mean.  Obviously,  any  errors  in  the  algorithm  are 
independent  of  the  day-to-day  variability  of  the 
ionosphere,  and  thus  these  two  errors  can  he  added  in 
root  of  the  sum  of  the  square  errors  fashion.  T  he  overall 
algorithm  improvement  obtained  in  this  manner,  from  the 
■191)  station  months  of  data,  was  approximately  55 
percent.  Thus,  these  initial  algorithm  tests  showed  that 
our  goal  of  a  50  percent  rms  ionospheric  correction  has 
been  met. 


ALGORITHM  LIMITATIONS 

T  he  algorithm  described  here  is  a  severely  truncated 
version  of  a  much  larger  empirical  model  of  ITX' 
developed  by  llent  |9|.  over  17  years  ago.  While 
significantly  more  actual  THC  data,  along  with  additional 
bottomsidc  and  topside  ionospheric  profile  measurements, 
are  available  than  at  the  time  of  the  development  of  the 
bent  model,  no  attempts  have  been  made  to  improve  the 
present  state  of  world-wide  empirical  TLX'  models.  T  his 
is  because  the  use  of  even  a  state  of  the  art  ionospheric 
model  would  not  get  around  the  approximate  ?.()  to  25 
percent  day-to-day  variability  of  the  ionospheric  THC 
which  the  single-frequency  user  must  live  with. 

T  he  high  latitude  regions  deserve  some  additional 
mention,  firstly,  the  auroral  regions  are  characterized  by 
high  variability,  with  little  chance  for  adequately 
modeling  departures  in  THC  from  monthly  average 
behavior,  except  in  a  statistical  sense.  Secondly,  most 
active  auroral  energy  inputs  occur  during  nighttime 
periods  when  the  background  absolute  TLX'  values  are 
considerably  lower  than  during  the  daytime  periods  when 
THC  is  higher.  Auroral  disturbances  nave  been  seen 
which  have  produced  ITX'  values  which  rival,  and 
sometimes,  even  exceed  midlatitude  midday  values,  but 
these  are  likely  infrequent  in  occurrence,  only  predictable 
in  a  statistical  sense,  and  cannot  be  used  to  improve 
single-frequency  user  ionospheric  range  errors. 
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CONCLUSIONS 

Tliis  algorithm  has  boon  implemented  in  a  Mightlx 
different  form  lor  use  with  the  Til’S  satellites.  Details  of 
its  implementation  in  a  single-frequency  user  set  are 
desaihod  in  detail  in  the  NAVST  AR  GPS  "Interface 
Control  Ooaiment''  |l()|.  It  is  relatively  easy  to  obtain  a 
SO  |H‘iveni  rms  ionospheric  error  reduetion  with  this 
algorithm  The  triek  is  in  the  single-frequency  user  being 
able  to  withstand  those  infrequent  times  when  a  large 
deviation  from  average  ionospherie  behavior  will  he 
eneountered.  If  a  single-frequency  user  cannot  tolerate 
those  .are,  but  real,  times"  when  tlfe  ionosphere  differs  by 
substantial  amounts  from  its  average  behavior,  then  that 
user  must  seriously  consider  opting  for  a  dual-frequency 
(JI’S  receiver  which  has  been  demonstrated  to  remove 
ionospheric  errors  to  a  level  comparable  with  other 
sources  of  (il’S  system  errors  1 1 1 1. 
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